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V 

Spreading  modulation  waveforms  are  discussed  In  general  terms,  and  a 
classification  system  Is  described.  Several  well-known  examples  are  given, 
and  a  new,  standard  set  of  spreading  modulations  Is  Introduced.  Spectra  of 
all  the  waveforms  are  given,  and  spectral  properties  are  discussed  In  relation 
to  system  performance.  It  Is  shown  that  these  properties  are  essentially 
determined  by  three  waveform  parameters,  for  a  large  class  of  spreading 
modulations. 
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1.  INTRODUCTION 

A  random  spreading  modulation  is,  basically,  a  means  of  complicating  the 
signal  of  a  communications  system  in  order  to  make  that  system  more  difficult 
to  jam.  Like  any  other  electronic  counter-counter  measure  tactic,  it 
is  a  game,  whose  object  is  to  maximize  the  cost  to  the  jammer  of  achieving  a 
given  degradation  to  one's  system.  Absolute  immunity  to  jamming  generally 
cannot  be  achieved,  except  possibly  temporarily,  until  the  opponent  figures 
out  what  you  are  doing.  The  tactic  of  spreading  is  approached  here  from  the 
point  of  view  that  the  adversary  knows  exactly  what  you  are  doing,  and  is  free 
to  do  his  best  to  degrade  your  performance. 

The  complication  introduced  by  random  spreading  has  two  essential 
features:  a  randomness  in  the  waveform  modulation  which  greatly  reduces  the 
effectiveness  of  repeat  jamming,  and  a  spreading  of  the  instantaneous 
spectrum,  which  increases  the  available  processing  gain  against  a  given  level 
of  total  jamming  power.  Simple  frequency  hopping,  which  increases  the  system 
bandwidth  without  changing  the  width  of  the  instantaneous  spectrum,  is  a 
related,  and  complementary  technique,  which  is  not  discussed  further  in  this 
study, 

A  simple  example  will  illustrate  the  application  of  a  random  spreading 
modulation.  Suppose  a  communications  link,  designed  for  a  friendly 
environment,  makes  use  of  a  signal  which  consists  entirely  of  simple  pulses, 
organized  in  a  way  which  permits  both  synchronization  and  data  transmission. 
The  pulses  all  have  the  same  shape  (say  rectangular)  and  the  same  carrier 
frequency,  and  occur  in  a  burst  of  some  fixed  pattern  for  synchronization. 

Data  may  be  transmitted  by  a  systematic  alteration  of  some  pulse  parameter, 
such  as  its  transmission  time,  as  in  pulse-position  modulation.  To  add  some 
anti-jam  capability  to  such  a  system,  the  simple  pulse  may  be  divided  into  N 
equal  sub-pulses,  or  chips,  with  the  carrier  phase  in  each  chip  either 
reversed  (i.e.,  changed  by  180°)  or  unaltered  in  accordance  with  an  N-bit 
pseudorandom  binary  sequence.  This  is  binary  phase-shift  keying,  BPSK,  one  of 
the  simplest  of  the  spreading  modulation  techniques.  By  the  use  of  BPSK,  the 
pulse  bandwidth  is  increased  by  a  factor  of  N,  which  is  one  of  the  objectives 
of  the  spreading  technique. 
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When  BPSK  Is  used  In  this  way,  the  pseudorandom  bit  sequence  must  be 
known  to  the  receiver,  where  a  filter  matched  to  the  spread  pulse  Is 
Implemented.  It  Is  a  major  feature  of  systems  utilizing  random  spreading 
modulations  that  the  spreading  bit-sequence  Is  known  to  the  receiver,  although 
any  modulation  technique  normally  used  for  the  transmission  of  Information 
bits  could.  In  principle,  be  used  as  a  spreading  modulation.  Thus  random 
spreading  Implies  a  different  receiver  structure  than  data  demodulation  of  the 
same  waveform,  and,  most  Important,  a  different  criterion  for  performance. 

It  Is  not  Intended  that  the  bit  sequence,  or  code,  used  In  one  pulse  be 
used  for  all,  but  rather  that  this  code  be  continually  changed  from  pulse  to 
pulse.  In  effect,  both  sender  and  receiver  must  have  the  ability  to  generate 
the  same,  unlimited  sequence  of  bits  In  time  synchronism,  successive  groups  of 
N  bits  being  used  to  “code"  successive  signal  pulses.  Of  course,  th?  sender 
must  be  able  to  modulate  his  pulses  with  these  ever-changing  code  sequences, 
and  more  significantly,  the  receiver  must  be  able  to  configure  a  filter, 
matched  according  to  some  criterion,  to  each  expected  coded  pulse.  For  the 
spreading  modulations  studied  here,  this  last  requirement  presents 
difficulties  of  widely  varying  degree.  Such  Implementation  Issues  will  have  a 
decisive  effect  as  the  choice  of  modulation  In  practice,  but  hardware 
techniques  change  rapidly,  and  no  attempt  Is  made  here  to  assess  the  value  of 
a  modulation  scheme  from  a  hardware  point  of  view. 

The  example  just  given  Is  perhaps  the  simplest  scheme  that  has  the 
requisite  properties  of  a  random  spreading  modulation.  A  very  general 
definition  would  consider  any  scheme  which  mapped  a  finite  sequence  of  bits 
Into  a  finite  segment  of  waveform  to  be  a  possible  spreading  modulation 
technique.  We  restrict  this  definition  by  the  application  of  two  constraints, 
one  of  practical  Importance,  and  one  of  a  simplifying  nature.  In  order  to 
define  a  workable  class  of  waveforms  for  analysis. 

The  first  constraint  Is  to  limit  ourselves  to  constant -envelope 
waveforms.  The  major  reasons  for  this  requirement  are  practical,  allowing  the 
use  of  simple  power  amplifiers  In  transmitters  and  repeaters.  Constant 
envelope  also  allows  the  maximum  signal  energy.  In  a  given  signal  duration.  In 
the  face  of  a  peak  power  limitation.  These  reasons  are  often  compelling  In 


real  system  design,  and  many  effective  spreading  modulations  have  been  devised 
within  the  constant-envelope  constraint.  Other  waveforms  may  be  converted  to 
constant-envelope  form  by  clipping,  either  deliberate  or  inadvertent,  but  the 
original  waveform's  spectral  properties  are  usually  not  improved  in  this  way. 
Subsequent  filtering  will  restore  some  amplitude  modulation,  hence  it  is 
desirable  to  study  constant-envelope  modulations  which  meet  the  spectral 
occupancy  limitations  within  which  the  system  must  be  designed  in  the  first 
place.  Spectral  occupancy  restrictions  play  a  key  role  in  spreading 
modulation  performance  evaluation,  and  the  general  goal  is  to  maximize  anti-jam 
processing  gain  while  operating  within  a  given  bandwidth  or  spectral  window. 

Since  the  waveform  is  to  have  constant  envelope,  it  can  be  described  by 
its  phase  or  frequency  modulation.  Our  second  constraint  requires  that  the 
waveform  be  segmented  into  equal  chips,  that  the  phase  or  frequency  modulation 
imposed  within  a  given  chip  be  one  of  a  small  number  of  possibilities,  and 
that  the  choices  made  in  the  successive  chips  depend  in  a  simple  way  on  a 
pseudo-random  bit  sequence  of  a  non-repeating  nature. 

We  tend  to  think  of  the  waveform  as  built  of  pulses,  as  in  our  simple 
example,  each  pulse  requiring  N  input  bits  for  its  specification.  But,  in  the 
analysis  of  spectral  properties,  we  treat  N  as  large,  and  it  is  not  important 
that  the  waveform  actually  be  segmented  into  separate  pulses,  hence  our 
results  will  be  valid  for  continuous  waveforms  as  well.  Our  second  constraint 
will  be  made  more  specific  in  Section  II,  where  a  classification  system  is 
introduced  and  the  method  of  spectral  analysis  described. 

If  a  signal  structure  is  indeed  built  out  of  pulses,  each  of  which  is 
modulated  by  a  random  spreading  modulation,  then  other  anti-jam  techniques  can 
still  be  used  as  well.  For  example,  pulse-to-pulse  frequency  hopping  can  widen 
the  system  bandwidth  still  further,  and  randomization  of  the  timing  of  pulses 
according  to  a  scheme  available  to  the  receiver  can  also  be  employed. 

Quantitative  measures  of  spreading  performance  have  been  discussed 
elsewhere^^^ ,  in  terms  of  efficiency  parameters  which  measure  the  loss  in 
processing  gain  (from  the  ideal  value  of  time-bandwidth  product)  which  results 
from  operating  in  a  fixed  band  against  an  optimizing  noise  jammer.  These 
parameters  depend  only  on  the  average  power  spectrum  of  the  modulation,  and 


they  favor  spectra  which  are  flat  within  an  allotted  band  and  fall  off  rapidly 
outside  that  band.  A  brief  derivation  of  these  performance  measures  is  given 
in  Appendix  A.  Interference  with  other  systems,  or  other  users  of  the  same 
system,  will  place  more  specific  requirements  on  the  spectral  density  of  the 
spreading  modulation,  in  terms  of  sldelobe  levels  and  spectral  decay  beyond 
the  nominal  band.  These  are,  basically,  spectrum  allocation  constraints,  and 
must  be  considered  separately  In  each  case. 

In  this  study  we  present  the  spectral  properties  of  various  groups  of 
spreading  modulation  waveforms,  without  judgment  of  their  relative  merits.  The 
actual  choice  In  a  particular  system  design  will  be  dictated  by  many  factors. 
Including  hardware  Issues  and  postulated  jamming  scenarios;  but  the  spectral 
properties  of  the  waveform  chosen  will  play  a  central  role  In  this  choice. 


II.  WAVEFORM  CLASSES  AND  METHOD  OF  COMPUTATION  OF  SPECTRA 

An  Initial  decomposition  of  spreading  modulations  Is  based  upon  two 
features:  the  constancy  (or  lack  of  It)  of  the  phase  within  each  chip,  and 

the  continuity  (or  lack  of  It)  of  phase  at  the  chip  boundaries.  Our  first 
class  consists  of  the  "constant  phase”  waveforms.  In  which  the  carrier  phase 
remains  constant  during  each  chip,  with  discontinuous  change  at  chip 
boundaries  permitted.  The  value  of  the  phase  In  a  given  chip  Is  determined  by 
a  particular  subset  of  the  code  bit  sequence  which  corresponds  to  the  pulse  In 
question.  In  the  simplest  case,  we  assign  the  bits  to  the  chips  In  a 
one-for-one  manner,  and  the  only  other  case  discussed  here  utilizes  two  bits 
per  chip  In  some  way.  Tlie  bit,  or  bits,  assigned  to  a  given  chip  can 
determine  the  phase  absolutely,  or  control  an  Increment  to  be  applied  to  the 
phase  of  the  preceedlng  chip.  When  two  bits  are  used  per  chip,  they  can  be 
two  "new"  bits  for  every  chip,  or  a  two-bit  window  can  slide  along  the  bit 
sequence,  controlling  phase  In  some  way.  These  possibilities  are  discussed  In 
detail  In  Section  III. 

The  other  major  class  of  modulations  combines  varying  phase  (throughout 

the  chip)  with  phase  continuity  at  chip  boundaries.  The  waveforms  are  most 

naturally  described  In  terms  of  the  possible  frequency  modulation  patterns 

employed.  One  can  have  a  repetolre  of^'^  possible  patterns  for  any  chip, 

selected  according  to  the  values  of  a  set  of  n  consecutive  bits  In  the  code 

sequence.  One  could  then  take  n  new  bits  for  each  chip,  or  slide  an  n-blt 

window  along  the  sequence,  yielding  a  correlation  between  successive  frequency 

modulation  patterns.  An  example  of  this  latter  group  called  "tamed  FM"  has 

(2) 

been  discussed  In  the  literature  as  an  Information  modulation,  but  the 
remainder  of  this  study  Is  devoted  to  the  so-called  "binary  FM"  waveforms,  of 
which  there  are  several  well-known  examples. 

The  binary  FM  waveform  uses  just  one  bit  for  each  chip,  hence  one  of  two 
frequency  modulation  patterns  Is  applied.  We  further  specialize  this  class  by 
the  assumption  that  one  of  these  patterns  Is  the  negative  of  the  other  (with 

respect  to  a  suitable  carrier  frequency),  so  that  the  frequency  modulation  In 

th  *  * 

the  n  chip  will  be  ±  <f((t-nA),  where  (|)(t)  Is  the  basic  frequency  modulation 
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pattern  of  the  waveform,  A  Is  the  chip  duration,  and  the  algebraic  sign  Is 
fixed  by  the  n^^  bit  of  the  spreading  code.  If  ^(t)  Is  a  constant,  the 
frequency-shift  keying  (FSK)  waveform  results,  the  simplest  example  of  a 
binary  FM  waveform.  Binary  FM  waveforms  are  discussed  In  detail  In  Sections 
IV,  V  and  VI. 

The  two  other  classes  allowed  by  our  original  division  are  excluded  from 
this  study  for  the  following  reasons.  The  first  class  would  consist  of 
constant-phase  waveforms  with  phase  continuity  at  chip  boundaries,  l.e., 
unmodulated  pulses,  and  these  are  of  no  Interest  as  spreading  modulations.  The 
remaining  class  allows  phase  discontinuities  at  chip  boundaries  along  with 
phase  variation  within  the  chip.  Ihe  phase  jumps  of  such  waveforms  Impart  to 
them  the  same  slow  decay  of  spectral  density  with  frequency  as  Is  exhibited  by 
BPSK  (namely  Inverse  square),  which  defeats  the  primary  purpose  of  the 
variable  phase  waveforms.  It  will  be  shown  In  Section  IV  that  the  asymptotic 
spectral  properties  are  determined  by  the  number  of  continuous  derivatives  of 
phase  at  chip  boundaries,  and  the  frequency  modulation  patterns,  4>(t),  of 
different  binary  FM  waveforms  are  chosen  with  this  property  In  mind.  To  allow 
phase  discontinuities  '■o  occur  (deliberately  or  accidentally)  would  ruin  the 
good  features  otherwise  attainable  In  the  waveform  design,  and  this  can  be 
clearly  demonstrated  In  specific  cases.  In  the  binary  FM  waveforms,  the 
frequency  modulation  function  Is  bounded,  and  the  resulting  phase  variation  Is 
continuous  over  the  entire  waveform.  If,  In  addition,  the  frequency 
modulation  function  Is  zero  at  the  beginning  and  end  of  a  chip,  then  the 
overall  frequency  modulation  of  the  waveform  Is  continuous,  regardless  of  the 
original  bit  sequence,  and  this  leads  to  faster  spectral  decay  with  frequency. 
Further  details  are  given  In  Section  IV. 

Our  basic  method  of  computing  spectra  begins  with  an  explicit  expression 
for  the  waveform  as  a  function  of  time  and  the  code  sequence.  Next,  we 
compute  the  Fourier  transform  of  this  time  function,  still  as  a  function  of 
the  code  bits.  The  square  of  the  Fourier  transform  gives  the  energy  spectrum, 
and  the  desired  power  spectrum  of  the  original  waveform  Is  obtained  by 
dividing  NA,  the  total  duration  of  the  signal,  which  Is  N  chips  long. 

Finally,  the  spreading  modulation  spectrum  Is  found  by  taking  an  ensemble 


average  over  the  bits  of  the  spreading  code.  These  bits  are  treated  as 
"purely  random",  i.e.,  independent  and  equally  likely  to  assume  either  value. 

The  "purely  random"  assumption  places  some  limitation  on  the 
applicability  of  our  results,  but  in  many  applications  one  strives  for  an 
approximation  to  this  quality  in  real  code  sequences.  This  is  not  a  basic 
limitation  of  the  method,  however,  which  would  still  be  applied  if  the 
statistical  properties  of  the  code  sequence  in  use  were  well  enough  known. 

The  spectra  of  individual  waveforms  will  vary  considerably,  depending  on 
the  actual  bits  of  the  code,  but  the  ensemble  averages  computed  here  are  still 
relevant  for  the  behavior  of  the  system,  so  long  as  code  sequences  change 
continually,  as  we  have  assumed  from  the  start. 

Other  methods  of  spectral  computation  are  equivalent,  but  tend  to  be  more 
difficult  to  apply.  For  example,  a  spreading  waveform  can  be  modelled  as  a 
stationary  random  process  by  extending  it  indefinitely  in  time  and  treating 
the  start  time  of  some  reference  chip  as  a  random  variable,  uniformly 
distributed  over  an  Interval  of  length  A,  while  the  chip  spacing  remains 
rigid.  One  can  then  compute  a  covariance  function  and  finally,  by  Fourier 
transformation,  a  power  spectrum.  Other  authors  have  used  Markov  chain  models 
to  describe  the  sequence  of  phase  variations  in  the  waveform,  and  these  too 
tend  to  become  complicated. 

The  essential  feature  of  the  method  used  here  is  a  representation 
theorem,  which  permits  the  expression  of  any  binary  FM  waveform  in  a  form 
formally  very  similar  to  the  representation  of  a  constant-phase  waveform.  This 
form  permits  easy  application  of  the  approach  outlined  above.  The 
representation  theorem  is  proved  in  Section  IV,  and  the  basic  spectrum 
calculation  is  given  in  Section  III. 
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III.  CONSTANT-PHASE  WAVEFORMS 

A  constant -phase  waveform,  N  chips  long,  Is  uniquely  specified  by  the 

sequence  of  carrier  phases  0.,  6,,  ...,  6„  This  sequence,  in  turn,  is 

0  1  N“i 

determined  from  an  input  bit  sequence  according  to  an  algorithm  characteristic 
of  the  modulation  scheme.  The  modulation  waveform  Itself  can  be  written  as  a 
summation  of  adjacent  rectangular  pulses,  each  one  chip  in  duration,  and 
having  the  appropriate  values  of  phase.  The  rectangular  pulse  function  Is 
defined  as 

fl  ;  0  <  t  <  A 

0  ;  otherwise, 

and  the  complex  modulation  function  which  describes  this  waveform  Is  then 

N-1  10 

Z(t)  -  I  e  "  P  (t  -  nA)  . 
n-0 

This  function  vanishes  outside  the  Interval  0  <  t  <  NA,  and  the  underlying 
carrier  la  arbitrary  (although  practical  generation  techniques  often  make  use 
of  an  Integral  relationship  between  chip  duration  and  carrier  period). 

The  familiar  Fourier  transform  of  the  rectangular  pulse  is  denoted  ko(u), 
as  follows: 

=  i 

o 

-la)A/2  sln(MA/2)  . 

®  .  (ajA/2) 

Fourier  transform  of  the  modulation  Itself  Is 

NA  .  ^ 

/  Z(t)  e"^“*^  dt 
o 

N-1  10  nA+A  ,  ^ 

-  J]  e  “  /  P  (t  -  nA)  e”  “  dt 

n-0  nA 

N-1  1(0  -  nuA) 

A  k^(u)  I  e  " 
n-0 


In  terms  of  kQ(u),  the 
K(a))  - 
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The  quantity  |K(a))|2  is  the  energy  spectrum  of  the  waveform,  and 
(according  to  Parseval's  theorem) 

^  /  lK(ai)rd6i  -  /  lZ(t)rdt  -  NA  , 

Q 

2 

since  |Z(t)|  is  identically  unity.  Thus  |K((ij)|  /NA  is  the  desired  power 
spectrum,  whose  ensemble  average  is  given  by 

G(f)  -  Ik(w)I^ 


-  A|k^(<o)|2 


N 


N-1  i(e  -  nojA)  , 
lie"  |2 

n-0 


We  write  spectral  densities  as  functions  of  ordinary  frequency,  with  the 
understanding  that  oi  -  2irf.  The  overbar  denotes  ensemble  average,  which  can 
be  evaluated  after  the  specific  dependence  of  the  phase  sequence  on  the  code 
bit  sequence  is  given. 

2 

The  separation  of  G(f)  into  a  "pulse  factor",  in  this  case,  |ko(<*>)l  » 

a  code  factor  (the  ensemble  average)  is  typical  of  all  the  spreading 

modulations  considered  in  this  study.  The  code  factor  is  periodic  in 

frequency,  with  period  equal  to  the  chipping  rate,  “  1/A,  hence  the 

spectral  behavior  for  all  constant-phase  waveforms,  at  large  frequency  (l.e., 

for  from  band  center),  is  given  by  the  pulse  factor,  which  in  this  case  decays 
-2 

as  f  .  Because  of  our  normalization  of  the  modulation  waveform,  the  power 
spectral  density  Integrates  to  unity: 

7  0(f)  f  .  1. 


The  code  factor  is  expanded  as  follows: 


N-1 

n-0 


^l(0n-nu)A)|2 


1 


N-1 

V 


A— (N-1) 


where,  for  A  >  0, 


C  =  1  y  ^ 

A  ■  N  ® 

n-0 


When  A  <  0, 


1  N-1 

“  n-L 


,^(®n+A  “  ®n) 


and  It  follows  that 


and  that 


C.A  -  C,  . 


C  -  1. 

o 


The  numbers,  C  ,  form  the  truncated  autocorrelation  sequence  of  the  finite, 
^  10 

random  sequence,  e 

If  the  phases  are  statistically  Independent,  and  If 


for  C^,  and  the  spectral  density  will  be 


_  ^  sln^(a)d/2) 

(o)A/2)^ 

This  Is  the  case  for  BPSK,  since  each  phase  angle  Is  equally  likely  to  assume 
either  of  two  values,  180  degrees  apart.  If  the  code  bit  sequence  Is  written 
b^,  bj^ . each  b^  Is  a  binary  variable  taking  only  the  values 

±1,  then  the  BPSK  phase  sequence  can  be  defined  by  the  statement 


-  0, 


then  all  the  C.  will  vanish,  except 
given  by  the  pulse  factor  alone: 


G(f)  -  Aik  (u)!' 
o 


for  0  <  n  <  N-1.  The  same  power  spectral  density  describes  the  QPSK  waveform 

(quadrlphase  shift  keying),  at  the  same  chipping  rate,  since  QPSK  can  be 

described  In  terms  of  a  sequence  of  2N  bits  as  follows: 

10  , 

'  <'’2n 


Four  phase  values  are  possible,  but  the  average  value  of  e^®^  is  still  zero. 

If  the  phases  are  independent,  but  if  e^®”  is  not  zero,  for  example  if 

e  *■  P  > 

then  lines  appear  in  the  spectrum.  It  is  not  hard  to  show  that,  for  large  N, 

2 

the  fraction  1  -  |p|  of  the  power  is  found  in  a  spectral  density  just  like 

2 

that  of  BPSK,  while  the  remaining  fraction,  |p|  ,  is  in  a  line  at  to  "  0  (band 
center).  This  kind  of  spectral  density  is  undesirable  in  a  spreading 
modulation,  hence  this  case  is  not  discussed  in  further  detail. 

A  number  of  constant -phase  waveforms  exhibit  correlation,  at  least  over 
adjacent  chips.  An  example  is  SQPSK  (staggered,  or  offset  QPSK),  in  which  the 
correspondence  of  code  bits  and  waveform  phases  is  best  shown  by  a  diagram; 


-2  ®n-l 


^2 


®n+l  ®n+2 


Here,  I  and  Q  refer  to  in-phase  and  quadrature  signals,  which  are  separately 
modulated  by  alternate  bits,  as  shown,  and  then  combined  to  form  the  resultant 
phase  sequence.  From  this  diagram,  we  can  write 

e  n-1  .  (b  +  lb  )  , 

✓2  "  "  ^ 
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e  ^  w 


‘"n  +  “'n.l)  • 


Only  adjacent  phases  are  correlated,  and  we  have 


,l(ea-an-l)  .  I  (b„  +  lVi)(t.„  - 


since  bn  bn+i^  -  bn  bn^^  ■  0,  for  It  +  0,  and  b^  ■  1.  Similarly  we  have 


,i(®n+l  “  ®n)  _  i. 


Cb  +  lb  ..)(b_  -  lb 


2  '  n+2  ’  ‘  n+l'^  n  "“n-l'  2  ’ 

and  this  equation  Is  thus  true  for  both  even  and  odd  values  of  n.  We  have 
therefore  found  that 

-  0  ,  Ul  >  1 


C  -  1  y  -  1  Ihi  _  c 

^1  N  ®  2  N  -1  • 

n«0 

For  large  N,  which  Is  always  the  case  of  Interest  here,  we  can  take  C+j  >  1/2 
and  thus 

G(£)  -  Aik^(a))i2  {1  +i  +  e^“^} 


.  sin  (uA/2)  .  _  .X 

A  - ^ ^  (1  +  cos  uA) 

(ajA/2)^ 


sln^((aA/2)  co8^(a)A/2) 
(uA/2)^ 


G(f)  -  2A  . 

(uA)^ 


In  this  case,  the  result  Is  the  same  as  the  spectrum  of  BFSK  or  QPSK  at  one- 
half  the  chipping  rate.  In  terms  of  our  diagram,  It  means  that  the  Q-channel 
signal  could  be  advanced  (or  delayed)  by  one  chip  length,  thus  aligning  the 
bits  to  produce  conventional  QPSK  with  chips  twice  as  long  as  the  original 
ones,  all  without  changing  the  spectrum. 
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A  second  example  exhibiting  one-chip  correlation  Is  Illustrated  In  the 
following  diagram: 

Q  t)n-3  bn_2  bn_i  b^ 

I  bii_2  bji-i  bji  bQ^.j 

®n-2  ®n-l  ®n  ®n+l 

We  call  It  PQPSK  ("Poor-man's  QPSK"),  and,  like  SQPSK,  the  bit  rates  and  chip 
rates  are  equal.  From  the  diagram, 

e  n  .  (b  +  ib  ) 

✓2  " 


'  -I  ‘Vi  +  -  “n-I*  ■  *'2  . 

For  large  N, 

Cl  -  1/2  ,  C_i  -  -  1/2, 

and  the  spectral  density  Is 

G(f)  -  A|k^(a))|^  {1+1  (e‘^“^  -  e^“^)} 

If  I  and  Q  channels  are  Interchanged  In  this  example,  the  second  factor  In  the 
spectral  density  Is  changed  to  (1  -  sin  uA). 

A  final  example,  which  has  been  called(^)  UPSK  (unidirectional  PSK) , 
shows  correlation  of  phase  over  two  chips.  The  defining  diagram  Is 


I  bn-l 


-2 

0n-l 

9n 

0n+l 

r 

0n+2 

1 

1 _ 

and  we  have 


9  1 

e  -  _L_  (b  ,  +  lb„  , ) 

n-1  n-1' 


1  1 

e  “-I  -  -i-  (b  +  ib  ) 
✓  2  “  "  ^ 


-  -!-  (b^  +  lb„) 

e  n  n 


e  ''+1  -  -L  (b  +  ib  ) , 
^2  ^■•■1  " 


etc.  We  evaluate 


^<Vl  -  ®n-2^ 


y  (1+1) 


1(0  -  0  ,) 
'  n  n-l"^ 


-^(1+1) 


-  I  (1  +  1), 


and  80  on,  and  also 


‘"’a  -  \-2>  „ 

e  -  0 


•  Vi>  ,  _ 


This  last  sequence  continues  to  alternate  between  the  values  0  and  1/2,  and  we 
find  (for  large  N); 


Cl  -  -i  (1  +  1)  ,  -  1/4  , 
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and  hence,  after  some  simplification, 

G(f)  -  2A  (1  +  sin  uA)  . 

(ojA)^ 

Again,  Interchange  of  I  and  Q  channels  in  the  definition  replaces  (1+sln  uA) 
by  the  factor  (1-sin  ojA). 

It  is  convenient  to  replace  frequency  by  an  angle  variable,  6,  defined  as 
follows : 

0  ■  ojA  2TrfA  -  Zirf/fchip* 

and  to  use  a  corresponding,  dimensionless  spectral  density: 
g(0)  =  G(f)/A. 

Then  g(0)  is  normalized  according  to 

CO 

^  /  g(0)  d0  -  1  . 


These  spectra  are  shown  In  Fig.  III-l,  where  spectral  density  is  plotted 
on  a  linear  scale.  The  plots  are  made  from  the  g(0)-formulas,  but  with  the 
abscissa  labelled  by  the  ratio  “  0/21:.  The  spectral  differences 

here  are  not  significant,  from  a  spreading  point  of  view,  and  the  disadvantage 
shared  by  all  Is  the  Inverse-square  decay  of  spectral  density  with  frequency. 
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UPSK 


Another  variation  of  the  constant  phase  theme  assigns  statistically 
Independent  Increments,  to  the  Initial  phase  of  each  chip  of  the  waveform 

Then 


and 


n 


1  -  9  +  y 

n  o 

» 

m 

m-1 

'•  n+1  n-* 

n+1 

-  exp(i  1  ^  ] 

1  . 

m-n+1 

If  the  Increments,  form  a  stationary.  Independent  sequence,  determined  by 
the  Input  bit  sequence,  then 


where 


> 


P  = 


The  code  factor  of  the  spectral  density  Is  then  given  by 


1  +  2Re 


N-1 

I 


1-1 


N-1 

N 


1  -ilaiA 
P  e  » 


which  becomes.  In  the  limit  of  large  N, 

_ 1  -  Ipl^ _ 

1  +  Ipl^  -  2Re(p  . 

In  general,  the  correlation  Introduced  by  specifying  increments  (Instead  of 
phases  directly)  leads  to  ripples  In  the  spectral  density,  as  given  by  this 
code  factor.  If,  however,  p«0,  then  the  BPSK  spectrum  reappears,  no  matter 
how  many  bits  are  used  to  specify  each  phase  increment. 

An  Interesting  case  occurs  if  each  phase  increment  is  either  ±\|), 
according  to  the  corresponding  code  bit.  We  can  write 


and  evaluate 


The  spectrum  of  any  waveform  featuring  Independent  phase  incremants  can  ba 
matched  by  one  of  these  special  waveforms,  with  appropriate  choice  of  ip.  The 
code  factor,  expressed  In  terms  of  ip.  Is 

sln^ip 

2 

1  +  cos  Ip  -  2cosip  cos((i}A) 

In  this  case,  and  we  shall  meet  this  factor  again.  In  a  similar  context.  In 
connection  with  binary  FM  waveforms. 

In  terms  of  0-aiA,  the  complete  spectral  density  for  this  waveform,  which 
we  call  IPSK  (Incremental  phase  shift  keying)  Is 

sln^ip 

2 

1  +  cos  Ip  -  2cosip  COS0 
Linear  plots  of  g(0)  for  a  range  of  ip-values  are  given  In  Fig.  III-2.  As  ^■►O, 
the  spectral  density  approaches  a  line  at  zero  frequency  (band  center),  and 
when  ip+ir,  lines  appear  at  frequencies  ±fchlp/2  on  either  side  of  band  center. 
The  relative  density,  G(f)/G(0)  Is  shown,  for  a  range  of  ip-values.  In  Fig. 
III-3,  expressed  In  dB.  The  null  positions  are  Independent  of  ip,  ana  t'*2  slow 
decay  with  frequency  Is  evident.  It  should  also  be  noted  that  sldelobe  levels 
Increase  with  increasing  values  of  ip  In  this  range. 


g(e)  - 


sin (0/2) 
072 
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Fig.  III-2.  IPSK  spectra 


lllT/16 


IV.  BINARY  FM  WAVEFORMS 


A  binary  FM  spreading  modulation  Is  completely  specified  by  Its 

characteristic  frequency  modulation  function,  <^(t),  which  Is  defined  only  over 
the  Interval  0<t<A,  one  chip  In  duration.  During  the  n^^  chip  of  the 
waveform,  which  starts  at  t-nA,  the  Instantaneous  frequency  Is  defined  to  be 

4)(t)  -  b^  ^(t-nA), 

where  bu  Is  the  n^^  bit  of  the  spreading  code.  The  Initial  phase  (at  t«0)  Is 
taken  to  be  zero,  so  that  the  modulation  waveform  Is 


where 


Z(t) 


l*(t) 


t  • 

<&(t)  -  /  <&(s)  ds. 

o 


The  last  chip  ends  at  t-NA,  and  Z(t)  Is  taken  to  be  zero  outside  the  Interval 
0<t<NA. 

The  characteristic  phase  variation  Is 

t  . 

•Kt)  -  /  41(8)  ds, 

o 

and  we  define 

4/  =  4(A). 


This  parameter,  1^  of  basic  Importance  In  the  study  of  these  waveforms;  It 
represents  the  magnitude  of  the  phase  change  which  takes  place  across  every 
chip.  The  actual  phase  Increments  accumulate,  with  appropriate  signs,  so  that 
the  waveform  phase  at  time  t-nA,  which  we  call  0^,  Is  given  by 

0  ■  *(nA)  ■  (b  +  b,  -f  ...  b  ,)4. 

n  '  '  o  1  n-1  ^ 

Of  course,  0^  «  0,  and  within  the  n^'^  chip,  we  have 


®(t)  “  0JJ  \  4(t-nA)  , 


where 


nA  <  t  <  (n-H)A 


Ve  assume  now  that  slnip  does  not  vanish,  deferring  discussion  of  this 
special  case  (where  ij;  Is  an  Integral  multiple  of  w)  until  later.  With  this 
restriction  we  Introduce  the  "pulse  function",  P(t),  which  Is  non-zero  only 
over  an  Interval  two  chips  long,  by  the  equations 

/  csct|»  slntij>-4>(t)] ;  0<t<d 

P(t)  -  ^  csc<»  sln[4)(t  +  A)];  -A<t<0 

\  0  ;  otherwise  . 

Note  that  P(t)  Is  continuous,  being  zero  at  t  -  ±A  and  unity  at  t-1,  but  not 
necessarily  symmetric.  Finally,  consider  the  waveform 

N  10 

Z  (t)  =  I  e  “  P(t  -  mA)  , 
m-0 

where  0  ■  (b  +  ...  +  b  ,  )ib  as  before, 

m  o  m-1  ^ 

For  a  time,  t,  within  the  n*^^  chip  (l.e..  If  nA  <  t  <  (n+l)A),  only  the 
terms  m»n  and  m»n+l  of  the  sum  defining  Zj(t)  contribute,  and  then 

10  10 

Zj(t)  -  e  "  P(t-nA)  +  e  P(t-nA-A) 

10 

e  ”  sln[i|»-^(t-nA)] 
sln[()i(t-nA)]}. 

+  lb  lb 
n^ 

> 

and 

lb  Ip 

e  ■  cos  (b^  Ip)  +  1  8ln(b^ip) 

-  cosip  +  Ib^  slnip, 

since  bn  Is  a  binary  variable.  We  therefore  obtain 


■  CSClp  { 


10 


+  e 


n+l 


But 


10  10 
n+l  n 

e  -  e 
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le 

Zj^(t)  -  csc\|/  e  ^  {sin(i(;-4i(t-nA)]  +  (co84M-ib^siin|))sin[<|)(t-nA)]}  , 
or 

19 

Zj(t)  -  e  ”  {co8[4i(t-nA)]  +  ib^  8in[4)(t-nA) } 

10  +  lb  (|.(t-nA) 

n  n  ^ 

■  e 

-  Z(t) 

Thl8  equality  of  Z(t)  and  Z^Ct)  holda  within  every  chip,  8lnce  the  Incluelon 
of  the  term  m>N  In  the  8um  defining  Z^(t) . valldatea  our  derivation  for  any 
value  of  n,  from  zero  through  N-1  (note  that  0^  Is  just  the  phase  of  Z(t)  at 
time  t«NA).  At  chip  boundaries,  Z(t)  and  Z^(t)  are  still  equal,  since 

10  16 
Zj(nA)  -  e  ”  P(o)  -  e  "  . 

The  modulation  function,  Zj(t),  Is  zero  when  t  <  -A  and  when  t>(N+l)A, 
and  thus  Z(t)  and  Zj(t)  differ  only  during  the  two  Intervals,  each  one  chip 
long,  which  precede  and  follow  the  original  time  Interval  over  which  Z(t)  Is 
non-zero.  The  new  function  Is  not  constant-envelope  In  these  "extra  chips,” 
but  for  large  N  their  presence  cannot  have  a  large  effect  on  the  spectral 
properties  of  Zj^(t),  and  henceforth  we  use  the  definition  of  Zj^(t)  as  a 
representation  of  Z(t),  dropping  the  subscript.  From  another  point  of  view, 
the  use  of  the  sum  for  Z(t)  Is  not  a  serious  approximation  because  Z(t),  as 
originally  defined.  Is  only  an  Idealization  of  the  waveform  likely  to  be 
produced  In  practice. 

For  many  purposes.  Including  the  computation  of  spectra,  the  sum 
representation  of  Z(t)  Is  very  much  more  convenient  than  the  original 
definition.  In  fact,  the  general  derivation  of  Section  II  needs  only  one 
change  to  apply  to  binary  FM  waveforms.  That  change  Is  to  replace  Po(t)  by 
P(t),  and  the  transform,  k^(oj),  by  the  new  transform 


Then 


-  ,  N  TTe  -  nuAj  ^ 

G(f)  -  a|1c(o»)|2  .  ^  I  ^  e  "  I 

n-0 

a  product  of  pulse  factor  and  code  factor,  as  before. 
The  code  factor  Is  simple,  since 


-  0  J  l(b  +  ...  +  b  ,  lip 

'■  n+l  n-'  '•  n  n+1-1-'^ 


-  e 


lb  i|) 

■  e 


lb 

e 


n+l-l 


m  (COS^)^  . 

This  follows  from  the  Independence  of  the  b^  and  the  Identity 


Ibtb  - 

e  -  cos^  +  lb  slnij;  ■  cosip  , 

for  any  binary  variable.  For  large  N,  we  get 

-  (cosip)^ 

for  1>0  and 

111 

S  -  C_^  -  (cos<>)'*' 

for  negative  1.  The  rest  of  the  evaluation  of  the  code  factor  has  already 
been  carried  out  In  Section  III  (waveforms  with  Independent  Increments),  and 
we  find 


1  N  1[0  -  nuAj  - 

^  I  I  e  ”  I' 

n-0 


N 

-  I 

1— N  * 


-iJtuA 


sln^ip 


1  +  COB  ii  -  2cos\p  cos 


Our  derivation  is  a  bit  heuristic,  but  the  result  is  easily  established  by  a 
careful  analysis,  so  long  as  |costp|<l. 

This  code  factor  is  unity  for  any  binary  FM  waveform  in  which  \|/  is  an  odd 
multiple  of  ir/2,  and  most  of  the  modulations  suggested  for  spreading  purposes 
share  the  value  t|;-ir/2.  It  turns  out  that  this  choice  leads  to  the  simplest 
receiver  design,  if  the  matched  filter  is  realized  passively,  but,  as  we  shall 
see,  other  values  of  tp  are  desirable  for  shaping  spectra  for  maximum  spreading 
effectiveness. 

The  asymptotic  behavior  of  the  spectral  densities  is  determined  by  the 

pulse  factor,  as  it  was  for  constant-phase  waveforms,  but  much  greater  variety 

is  now  possible.  The  discontinuous  pulse  function,  F^(t),  led  to  an  Inverse 

square  decay  In  the  former  case,  but  all  binary  FM  waveforms  have  spectral 

-4 

densities  which  fall  off  at  least  as  fast  as  f  (it  is  assumed  that  the 
characteristic  frequency  modulation  function,  ^(t),  is  bounded).  To  verify 
this  property,  we  write 

-i(i»t  1  d  -iojt 
®  “  (0  dt  ®  ’ 

and  carry  out  a  partial  Integration  in  the  defining  equation  for  k(u)): 

k(a))  -  I  /  P(t)  e'^^’^dt 
—A 

■hi 

—A 


(i>A 


/  P(t)  e~^“*^dt. 
-A 


The  Integrated  part  vanishes  because  F(t)  Is  continuous  throughout  the 
Interval,  Including  the  point  t-0  and  the  points  t-±A,  where  it  is  zero. 
However,  F(t)  may  be  discontinuous  at  these  points,  and  a  repetition  of  the 
procedure  yields 
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kCu)  -  -5-  /  P(t)  de  “ 

(D  A  -A 

- ^  /  P(t) 

0)  A  -A 

.  1  r^/ A  \  -IwA  *,  . , V  iuA 

+  -5-  {P(A-)e  -  P(-A+)e 

0)  A 

•  • 

-  P(0+)  +  P(0-)}  . 

•  • 

The  notation  P(t±)  stands  for  the  limit  of  P(t±e)  as  c*0,  and 

we  see  that  k((i))  varies  as  f“2  (hence  G(f)  goes  like  f“^)  If  P(t)  Is 
discontinuous  at  t-0  or  t-±0,  and  If  P(t)  Is  bounded  within  (-A,A). 

If  P(t)  Is  also  continuous,  the  argument  Is  repeated  again,  and  so  on,  so 
that  one  sees  that  the  asymptotic  properties  of  G(f)  are  directly  correlated 
with  the  smoothness  of  the  pulse  function,  P(t),  at  the  origin  and  at  the  ends 
of  Its  range. 

The  smoothness  of  P(t)  depends.  In  turn,  on  the  corresponding  properties 


of  <Kt).  For 

Instance,  In  the  Interval  -A<t<0, 

•  • 

P(t)  -  csc^  C08[(fr(t+A)]  ^(t+A)  , 

and  hence 

•  • 

P(“A+)  -  csc^  <>(04-) 

and 

•  • 

P(O-)  ■  cot(|)  (^(A-)  . 

Similarly, 

and 

•  • 

P(0+)  -  -cottp  ()»<04-) 

P(A-)  -  -csc<)  <fr(A-) 


If  P(t)  Is  to  be  continuous,  i.e.,  zero,  at  t«±A,  then  ^(t)  must  vanish  at 
both  ends  of  Its  range  of  definition.  This,  In  turn,  makes  P(t)  continuous  at 
t-0.  If  i^(t)  vanishes  at  t*0  and  t-A,  then 
•• '  •• 

P(-A+)  -  csci|)  (ji(0+) 

P(A-)  -  -csci|)  (|i(A-) 

and  so  on,  hence  the  controlling  parameter  In  this  matter  Is  the  order  of  the 
zero  of  (^(t)  at  t^O  and  t-A  (or  the  lower.  If  these  are  unequal). 

It  Is  easy  to  see  why  this  should  be  so,  by  considering  two  adjacent 

s 

chips.  At  the  chip  boundary  the  Instantaneous  frequency  changes  from  ±^(A-) 

•  s 

to  ±  ^(0+),  depending  on  the  relevant  bits,  hence  the  more  smoothly  ^(t) 
approaches  zero  at  the  ends  of  Its  range  the  smoother  will  be  the  transitions 
from  one  chip  to  the  next. 

To  summarize  this  relationship,  we  can  say  that  If  Is  the  lowest- 

order  phase  derivative  which  falls  to  vanish  at  t-0  and/or  t-A,  then  k(a))  will 
decay  as  f  ,  and  the  spectral  density  as  f”  ”,  at  great  distance  from 
band  center.  This  statement  Is  also  true  for  constant  phase  waveforms,  which 
correspond  to  n*0,  while  for  all  binary  FM  waveforms,  n>l.  The  parameter  n, 
along  with  ip,  has  a  major  effect  on  the  character  of  the  spectra  of  binary  FM 
waveforms. 

We  conclude  this  section  with  a  discussion  of  some  special  values  of  i(i, 
the  phase  shift  per  chip.  First,  let  iji-x/Z,  which  Is  typical  of  many 
waveforms  used  for  the  transmission  of  Information.  As  noted,  the  code  factor 
Is  unity  for  these  cases,  hence  the  spectral  density  of  the  waveform  Is 
Identical  to  that  of  the  pulse  factor  Itself.  The  basic  representation  also 
simplifies,  since 

.  lb 

for  any  binary  variable,  and  therefore,  for  n>0, 

10 

e  -1  bb,...b 

0  1  n-1 
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We  can  define  a  new  bit-sequence,  a^i,  as  follows: 
■  1- 


•l  ■  ‘o  ■  ‘o^o  • 

*2  ■  Vl  ■  ‘ih' 

a  *b  •••b  _*a  ,b 

n  o  n-1  n-1  n-1’ 

and  then  write 

N 

Z(t)  -  I  l”  a  P(t-nA)  . 
n-0 

The  a-sequence  Is  purely  random  If  the  b-sequence  Is,  and  the  b's  are 

recovered  from  the  a's  by  means  of  the  relation 

b  -  a  a  . 
n  n  n+1 

The  new  expression  for  Z(t)  Is  particularly  useful  for  the  study  of  techniques 
for  generating  these  waveforms  and  for  the  design  of  matched  filters  for  the 
corresponding  receivers.  The  case  -ir/2  leads  to  the  complex  conjugate  of 
this  representation  for  Z(t). 

When  tp-Kx,  the  whole  analysis  must  be  changed,  and  we  go  back  to  the 
original  definition  of  Z(t),  namely. 


Z(t)  -  e 


10^  +  Ib^  ♦(t-n&). 


for  values  of  t  In  the  n  chip.  But  now. 


lb<p  IbKx  „  /w  \  /  ,\K 

;  -  e  -  cos(Kx)  -  (-1)  , 


hence 


n  ,  ,.nK 
e  -  (-1)  , 


which  Is  Independent  of  the  bit  sequence.  Then 

Z(t)  -  (-1)"K  {cos(#(t-nA)l  +  Ibn  8ln(*(t-nA)] } 
within  the  n^^  chip,  and  the  I-component  Is  deterministic.  If  we  define 


A(t)  -  C08[<jl(t)]  PjjCt), 

B(t)  -  8ln[<|i(t)I  Pjj(t), 

where  P^Ct)  l8  the  rectangular  pulse  function  defined  In  Section  III,  then  we 
can  write 

N-1  „ 

Z(t)  -  I  (-1)”*^  A(t-nA) 
n-0 

N-1  nif 

+  1  I  b  B(t-nA) 

n-0 

=  X  (t)  +  lY(t)  . 

In  terms  of  the  pulse  function  Fourier  transforms: 
k^((o)  =  I  /  A(t)  dt 

-  j  /  co8t<^(t)]  e~^“^dt 

-A 

and 

«0 

ky(a))  =  I  /  B(t)  e~^“’'dt 

-  j  !  8lnI.Kt)]  e’^^^'dt  , 


we  obtain 


K((o)  -  /  Z(t)  dt 

N-1  „  4  * 

.1  /  \  r  /  ,\nK  -InujA 

-  Ak  (oj)  I  (-1)  e 

*  n-0 

+  Ak  ((0)  Y 
y  n-0 

=  K  (oj)  +  IK  ((o)  . 

A  y 

The  spectral  density  Is  nou 


0(f)  Ik(»)|2  -Ij  |K,(»)I^  ly»)|2  . 

In  the  first  term,  we  evaluate 

N“1  <nMA  ^  (wA4Kir)  sin  ^  (uA+Kir) 

I  (.i)"K  g-lna)A  _  ^  2  .  - 2 -  ^ 

n-0  sin  ((oA+Kir) 

and  hence 

,  sln^  (ojA+Ktt) 

-AlkCm)^  - ^ -  . 

N  sln‘‘  j  (wA-HCir) 

For  large  N,  the  second  factor  represents  a  sum  of  5-functlons  at  the  points 
ojA+Ktt  -  2LTr,  for  Integral  L,  and  G^(f)  becomes  a  line  spectrum.  In  the  usual 
case,  all  but  a  finite  number  of  these  lines  are  cancelled  by  zeroes  of  the 
factor,  k  (ui). 


The  Q-component  tero  represents  e  continuous  spectrua: 


-  4  lkM|2  1  lY  b 

^  ”  n-0 

-  A|ky(u))|^  , 

In  view  of  the  Independence  of  the  bits  In  the  code  sequence. 

Since  spectral  lines  are  very  undesirable  In  a  spreading  waveform,  the 
analysis  of  these  cases  Is  carried  no  further.  It  Is  worth  noting,  however, 
that  upon  squaring,  the  phase  modulation  of  a  waveform  Is  doubled,  so  that  a 
spreading  waveform  with  a  given  value  of  i()  will  look  like  a  waveform  of  the 
same  type,  but  with  double  the  ip-value,  after  passing  through  such  a  non¬ 
linearity.  Thus  a  waveform,  like  MSK,  with  will  yield  lines  In  the 

spectrum  If  squared,  and  this  could  be  undesirable  from  the  point  of  view  of 
signal  detection  and  frequency  determination  by  a  potential  jammer. 


V.  TYPICAL  BINARY  FM  WAVEFORMS 

The  simplest  binary  FM  waveform  Is  FSK  (frequency  shift  keying),  defined 
by  the  statement 

e 

<^(t)  -  constant. 

e 

The  constant  value  of  may  be  written  <|>/A,  and  the  carrier  frequency  Is 

Increased  or  decreased  by  this  amount  during  each  chip,  according  to  the  sign 

of  the  corresponding  code  bit.  Since  ^(t)  Is  never  zero,  we  know  that  the 

-4 

spectral  density  will  behave  like  f  (n-1)  as  f-M», 

The  FSK  pulse  function,  P(t),  Is  shown  In  Fig.  V-1  for  various  values  of 

<{).  The  special  case,  t(i-7r/2.  Is  the  familiar  MSK  (minimum  shift  keying) 

waveform,  distinguished  here  by  the  continuous  slope  of  P(t)  at  the  origin. 

All  binary  FM  waveforms  having  i|)“Tr/2  share  this  latter  property,  due  to  the 
factor  cott(i  In  the  expressions  for  P(0±).  However,  the  slope  discontinuity  of 
the  pulse  function  at  t-±A  Is  then  the  determining  factor  for  asymptotic 
spectral  behavior. 

The  Fourier  transform  Is  easily  evaluated.  After  noting  that 
P(t)  -  cscip  sln(i|;  -  , 

we  write 

A 

k(u)  _  .  2  I  8ln(i|»-  ip)cos  cot  dt 

o 

A 

-  j  {sln(c(i  -  4  +  cot)  +  sln(c|>  -  4  ’t’  ”  , 

o  A  Zi 


FSK 


2 

We  combine  |k.(a))|  with  the  code  factor  and  express  the  result  in  terms  of  the 
angle  variable,  &-a)A,  and  the  dimensionless  g(0)  -  G(f)/A,  as  before: 

/  -V  ,  .2  1  cos  e-cos 

8(0)  - 

{  9  -'I' 

Plots  of  this  function  on  a  linear  scale  are  given  in  Fig.  V-2  for  a  wide 
range  of  tj;  values.  The  effect  of  the  denominator  of  the  code  factor  is 
obvious,  as  it  was  in  the  case  of  IPSK,  and  as  \Ji  approaches  it  or  zero,  the 
spectral  density  acquires  lines. 

It  is  apparent  that  g( 0)  can  be  made  fairly  flat  in  its  "mainlobe",  by 
proper  choice  of  tji,  and  this  is  shown  in  Fig.  V-3,  where  g(0)/g(O)  is  plotted 
in  dB,  for  a  narrower  range  of  ilrvalues.  It  should  be  noted  that  the  first 
sldelobe  increases  steadily  with  increasing  \(i.  In  Appendix  A,  two  quantitative 
measures  of  "spreading  efficiency"  are  introduced,  and  it  is  shown  that  these 
quantities  are  optimized,  in  the  case  of  FSK,  for  \|r-values  near  those  which 
produce  good  "flatness"  of  the  spectral  density,  in  spite  of  the  effect  on  the 
first  sldelobe.  Those  performance  measures  deal  with  spectral  occupancy 
constraints  in  a  very  simple  way,  and  in  some  applications  the  sldelobe  level 
could  have  more  significance. 

The  expected  behavior  of  g(0)  with  large  0  is  shown  in  Fig.  V-4,  for  a 
wide  range  of  i|rvalues.  The  quantity  plotted  is  log  R(f),  where  R(f)  is  the 
fraction  of  the  total  power  which  lies  outside  a  band,  of  width  2f,  centered 
on  the  carrier.  In  other  words. 


where  0  ■  2xA.  These  curves  are  obtained  by  numerical  integration. 


1 


1+cos  -  2cos\J)  COS0 


3A 


chip 


spectra  of  FSK 


The  SFSK  waveform  (sinusoidal  FSK),  originally  Introduced(^)  In  the 
special  version  with  ^ir/2,  Is  defined  by  the  frequency  modulation  pattern 

♦(t)  -  -I  [1  -  coa(2irt/A)]  . 

•  •• 

Not  only  Is  <^(t)  zero  at  t-O  and  t-A  (unlike  FSK),  but  ^(t)  vanishes  at  these 

points  also.  Thus  n>3  and  the  spectral  density  will  vary  as  f  as  f-M*.  The 
SFSK  phase  variation  Is 

*(t)  -  I  -  -57  8ln(2itt/A)}  , 

which  has  the  property 

iKA-t)  -  ♦  -  i|i(t)  , 

as  a  consequence  of  the  symmetry  of  ^(t)  about  the  midpoint,  t»A/2.  From  the 
definition  of  P(t)  we  see  that  this  symmetry  of  +(t)  also  Implies  symmetry  of 
the  pulse  function  about  t«0.  These  symmetry  properties  are  shared  by  all  the 
binary  FM  waveforms  discussed  In  the. present  study. 

Since  F(t)  Is  even,  the  Fourier  transform,  k((i)).  Is  computed  from  the 
expression 

k(u)  .  I  sln[^  -  4(t)]  cos  (i)t  dt 

o 

■  "  (f  "  •l-*‘(2nt/A)] 

o 

+  sln[^  “  ^"a  ^  8in(2iTt/A)] }  dt. 

(4) 

Following  Amoroso  ,  the  Integrand  can  be  expanded  In  a  series  of  Bessel 
functions,  after  which  term-by-term  Integration  yields  a  rapidly  convergent 


series.  However,  k((D)  can  be  computed  equal  simply  by  numerical  Integration, 
since  the  range  of  Integration  Is  finite  and  the  Integrand  Is  oscillatory,  so 
that  high  accuracy  can  be  obtained  with  a  modest  number  of  points.  Numerical 
Integration  has  the  advantage  of  applying  equally  well  to  all  binary  FM 
waveforms,  and  the  procedure  Is  systematized  and  described  In  Section  VI. 

Spectra  obtained  In  this  way  are  shown  In  Figs.  V-5,  V-6  and  V-7,  which 
are  analogous  to  corresponding  plots  for  FSK.  The  effectiveness  of  \|;  In 
controlling  flatness  Is  again  seen  here,  as  well  as  the  anticipated  behavior 
at  large  values  of  f.  A  perhaps  unexpected  result  Is  the  relatively  large 
first  sldelobe,  which  Is  In  all  cases  greater  than  that  for  FSK  with  the  same 
i|^values.  This  appears  to  be  a  penalty  associated  with  Increasing  values  of 
the  Index  parameter  n,  as  we  shall  see  In  more  detail  In  Section  VI. 

Two  examples  of  waveforms  with  Index  parameter  n»2  (G(f)  ~  f”^)  are 
discussed  next.  We  call  them  waveforms  A2  and  B2,  and  they  have  been 
briefly  discussed  In  the  literature^ . 

Waveform  A2  Is  defined  by  the  frequency  modulation  pattern 

■{•(t)  -  8ln(Trt/A)  , 
with  the  corresponding  phase  variation 

•Kt)  -■!  [1  -  co8(Trt/A)]  . 

• 

Clearly,  n-2,  since  i^(t)  Is  Itself  the  highest  derivative  of  phase  which 
vanishes  at  both  ends  of  the  chip  Interval.  This  waveform  Is  similar  to  SFSK 
In  the  sense  that  Its  frequency  modulation  pattern  Is  expressed  In  terms  of 
trigonometric  functions,  but  It  Is  one  degree  less  smooth  at  the  chip 
boundaries.  Spectral  plots  on  linear  and  dB  scales  are  given  In  Figs.  V-8  and 
V-9,  respectively,  and  Fig.  V-10  shows  the  cumulative  spectra  In  the  form  of 
plots  of  log  R(f).  The  general  behavior  Is  similar  to  FSK  and  SFSK,  and  the 
magnitude  of  the  first  sldelobe  is  Intermediate  between  these  two. 

Our  second  example  having  n»2,  waveform  B^,  Is  defined  by  the  frequency 
modulation  pattern 
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Fig.  V-5.  SFSK  spectra 


6\li-  it/16 

Gff )  I  /'P*  13Tr/16 


Relative  spectra  of  SFSK 


i|^137r/16 


Fig.  V-9.  Relative  spectra  of  waveform 


TR-596(V-10) 


Fig.  V-10.  Cumulative  spectra  of  waveform  A, 


1(4i|>/A)(t/A)  ;  0  <  t  <  A/2 

(4t|;/A)  (l-  I)  ;  A/2  <  t  <  A  . 

This  function  has  a  triangular  shape.  Instead  of  the  half  sinusoid  of  waveform 
and  features  a  discontinuous  derivative  at  t-A/2.  This  discontinuity  is 
of  the  same  order  as  those  at  t-0  and  t-A,  and  therefore  does  not  alter  the 
expected  character  of  the  asymptotic  spectra.  The  phase  variation  is 

I2t|;(t/A)2  ;  0  <  t  <  A/2 

t|»[l-2(l-t/A)2]  ;  A/2  <  t  <  A  , 

and  the  spectra  are  presented  in  the  same  three  forms  as  before  in  Figs.  V-11, 
V-12  and  V-13.  These  latter  are  closely  comparable  to  the  plots  for  waveform 
^2* 

^  A  * 

The  normalized  frequency  modulation  patterns,  -^  ()>  ,  of  the  four  spreading 

modulations  discussed  in  this  section  are  plotted  together  in  Fig,  V-14,  and 

the  corresponding  pulse  functions  are  illustrated  in  Fig,  V-15  for  the  special 

case  t|)>ir/2.  In  this  latter  case,  we  have 

Icost((((|t|)]  ;  |t|  <  A 

0  ;  otherwise, 

which  is  direct  generalization  of  the  well-known  expression  for  MSK. 

The  four  waveforms  are  further  compared  in  the  Table  V-1. 

TABLE  V-1 

SIDELOBE  LEVELS 


Waveform 

n 

SL  (Tr/2) 

SL(t|;) 

FSK 

1 

-23.0  dB 

-18.2  dB 

5ir/8 

A2 

2 

-20.6  dB 

-17.6  dB 

19Tr/32 

B2 

2 

-20.0  dB 

-17.2  dB 

19Tr/32 

SFSK 

3 

-18.7  dB 

-16.3  dB 

37Tr/64 

The  column  SL(ir/2)  gives  the  magnitude  of  the  first  sidelobe  of  each  waveform 
for  the  choice  i|)-ir/2.  The  column  SL(<|;)  gives  that  sidelobe  level  for  a  value 
of  \|;  which  might  be  chosen  to  optimize  the  mainlobe  spectral  flatness,  and  the 
corresponding  value  of  ^  (essentially  an  eyeball  choice)  is  shown  in  the  last 
column.  First  sidelobe  levels  are  plotted  as  functions  of  in  Fig.  V-16. 


VI.  POLYNOMIAL  FM  WAVEFORMS 

It  Is  apparent  from  the  spectra  of  the  waveforms  studied  In  Section  V 
that  the  parameters  t|;  and  n  have  a  major  Influence  on  spectral  properties  and 
that  waveforms  such  as  A^  and  with  the  same  n-value,  have  very  similar 
spectra  for  equal  values  of  ip.  In  order  to  study  the  effects  of  these 
parameters  systematically,  a  special  set  of  waveforms  has  been  devised,  with 
one  member  for  each  value  of  n.  The  other  parameter,  i|>.  Is  simply  a  scale 
factor  for  the  characteristic  frequency  modulation  pattern,  as  usual. 

To  control  the  smoothness  of  (|>(t)  In  a  convenient  way,  we  make  use  of 
simple  polynomials,  and  postulate  the  general  form 


The  phase  function  Itself  will  be  proportional  to  t  near  t-0,  with  an 
analogous  behavior  near  t-A,  hence  the  first  phase  derivative  which  falls  to 
vanish  at  these  points  Is  the  n*”^.  Thus  the  subscript  correctly  corresponds 
to  the  value  of  the  Index  parameter  Itself.  The  frequency  modulation  patterns 
are  symmetric  about  t»A/2,  and  hence  the  pulse  functions  will  be  symmetric 
functions  of  t.  These  waveforms  are  not  the  same  as  the  polynomial 
modulations  discussed  by  Slmon^^\ 

The  phase  variations  are  obtained  by  Integration: 

.1.  ^  _  1  t. '  n“l  , 


(1 -f)  dt 


-  (1-x)""^  dx. 

o 


Since  'l'jj(A)  must  equal  ip,  the  constants  are  determined  by  the  requirement 


1  -  K  /  x""^  (l-x)“"^  dx 

n  ' 
o 


'.'.t 


1  2 

-Ik  f 

*>  n  J  '•A 


-1 


n-1 
)  du 


K  Or.  1 

-  /  (cose)^"^  de  . 


This  Is  a  knovm  Integral,  and  we  find 


K 


n 


(2n-l)l 

[(n-1)!]^ 


“  ( 


2n-l 

n-1 


)  . 


It  should  be  noted  that  for  n^l,  the  waveform  Is  simply  FSK. 

We  call  these  waveforms  "polynomial  FM”  waveforms,  and  denote  them  by  the 
symbol  PFM^.  The  normalized  frequency  modulation  patterns  for  the  first  five 
waveforms  are  compared  In  Fig.  VI-1*  In  Fig.  VI-2  we  compare  these  patterns 
for  PFM^,  Aj  and  while  Fig.  VI-3  compares  PFM^  with  SFSK.  From  the 
similarity  of  those  frequency  modulation  patterns  sharing  a  common  value  of  n, 
we  Infer  that  our  waveforms  will  provide  a  representative  set,  covering  the 
entire  class  of  binary  FM  spreading  modulations. 

The  computation  of  spectra  by  numerical  Integration  Is  facilitated  by  a 
change  of  variable  which  takes  advantage  of  the  symmetry  properties  of  the 
^(t).  For  any  frequency  modulation  pattern  symmetric  about  t-A/2,  we  will 
have  the  corresponding  property 

♦(A-t)  ~  i,  -  ♦(t), 

and  the  pulse  function  will  be  even  In  t.  From  this  last  fact,  we  can  write, 
as  before. 


or 


k(o)) 


2 

A 


A 

/  P(t)  cos  <i)t  dt  , 
o 


slnip  k(o)) 


sln[i|»-f(t)]  cos  ci)t  dt 
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o 


Now 


♦  ■»*)  -  I 

is  an  odd  function  of  x,  hence  we  Introduce  the  new  variable,  u,  by  aeans  of 
t  -  (1+u)  , 


This  variable  ranges  from  -1  to  +1,  and  we  can  put 
♦(t)  -  ♦(! +y^)  =  |tl  +  h(u)]  . 

The  new  function,  h(u),  is  an  odd  function  of  u  with  the  boundary  values 
h(0)  -  0  ,  h(l)  -  1, 

and  it  describes  the  phase  variation  over  half  a  chip  in  a  normalised  fora. 
In  terms  of  the  new  variable,  u: 

sin,  k((ii)  ■  /  (•Inf'l  *  f 

+  .in[|-|h(u)  -f-ifl])  J„  . 

We  now  separate  into  even  and  odd  functions  of  u: 
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.  .ln(i^)  CO.  [<*■(")  -  “'^“] 


/*+4JA-i  .  r\(lh(u)  -  UJAUi 

-  cos  (-5^)  sin  'a - ]  , 


.  rii  OlA  1./  X  to)AUi 
sln[^  ■  y  "  2  ■  ~J 


8ln  (-StliiA)  cos 


/-ib-aiAx  .  r'l'h(u)  +  oiAui 
-  cos  sin  - J 


The  odd  terms  vanish  upon  integration,  and  we  obtain  the  desired  result: 
sin4ik(a))  ■  sin(^‘**^)  /  cos  j  ju 


+  ,1„  (jeSi)  ;  CO.  (WU)  -  lOiOj  _ 


In  terms  of  6*uA  and  the  dimensionless  spectral  density,  g(9),  we  have 


g<0)  - 


[sln(i^)  f^(6)  +  .in  (^)  f_(e)] 
2 

1  +  cos  ij)  -  Icos^t'  cos 6 


where 


f  +  (0)  =,  f  COB  [- 


);h(u)  t  9ui 


These  last  integrals  are  easily  computed  numerically,  and  as  few  as  16  points 
yields  adequate  accuracy  for  most  purposes. 
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To  find  h(u)  for  a  given  waveform,  it  is  simplest  to  begin  with  its 
derivative,  h'(u),  which  satisfies 

i(t)  -  I  h'(u)  -  I  h'(u)  . 

For  the  PFM  waveforms,  we  write  h'n(u)  for  the  function  corresponding  to 
i^(t).  Since 

1  1  +  u 

A  "  2  ’ 

we  have 


n-1  ,  n-1 

Kq  (^)  (^)  -  h'nCu) 


K  <,  n-1 

h'  (u)  ■  — ~  (l  -  u^)  , 

n'  ,n-l  '' 

A 


The  corresponding  phase  function  is 


V“>  ?  Cl)  (-1)'  7 

a"  ^  i-O  *  o 


h  (u)  -  — u  y  (“"h  1"“  J 

n^  ^  ,n-l  .“ft  21+1 


It  can  be  shown  that  h^(l)-l,  as  required,  and  the  first  five  polynomials  are 
listed  here  for  reference: 


hj(u)  -  u 

h^Cu)  -  -I  (3  -  u^) 

hsCu)  -  I  (15  -  lOu^  +  3u^) 

h^(u)  -  (35  -  35u^  +  21u^  -  5u^) 

h5(u)  -  (315  -  420u^  +  378u^  -  180u^  +  35u®)  . 

For  comparison,  the  h-functions  corresponding  to  the  other  waveforms  studied 
In  Section  V  are  as  follows: 

Waveform  :  h(u)  ■  sin 

Waveform  :  h(u)  ■  u(2  -  lu|) 

SFSK  :  h(u)  -  u  +  —  sln(iru)  . 

w 

The  FFM  waveforms  exhibit  the  same  kind  of  spectral  behavior  that  we  have 
seen  In  FSK  and  the  other  examples  of  Section  V.  For  very  small  values  of  i(;, 
all  the  spectra  become  narrow  lines  at  band  center  (f-0),  and  as  ij;  approaches 
the  value  it,  In  all  cases  spectral  lines  appear  at  coA  -  ±ir  (l.e,,  f-± 

The  Interesting  region  lies  between  these  bounds  and  In  the  plotted  spectra, 

Is  confined  to  a  range  of  values  between  it/2  and  11  ir/16. 

The  plots  In  Figs.  VI-A  through  VI-13  show  spectra  for  waveforms  FFM^ 

(FSK)  through  FFM^,  The  FSK  spectra  are  Inculded  again  for  comparison.  For 
each  waveform,  linear  plots  out  to  relative  spectra  (In  dB)  out  to 

f-2f^jjj^p  are  given.  In  each  case,  \|>  ranges  In  Increments  of  it/32  from  the 
value  it/2  through  11  tt/16.  The  first  sldelobe  levels  are  shown,  as  functions 
of  In  Fig.  VI-14. 

The  concept  of  spreading  efficiency  Is  discussed  In  Appendix  A  as  a  measure 
of  the  processing  gain  attainable  against  an  optimizing  noise  Jammer,  compared 
to  the  nominal  tlme-bandwldth  product  available  to  the  system.  It  Is  shown 
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Fig.  VI-4.  PFM  (FSK)  spectra 


Relative  spectra  of  PFM.  (FSK) 


tra  of  PFM, 


Fig. 


there  that  FSK  waveforms  exhibit  a  maximum  In  spreading  efficiency  for  a 
particular  value  of  and  that  this  value  corresponds  closely  to  the 
"flattest”  spectral  density.  The  same  behavior  Is  shown  by  all  the  PFM 
waveforms,  and  for  each  of  these  a  \((-value  has  been  chosen  which  Identifies 
the  flattest  member  of  the  set  of  spectra,  as  Judged  by  eye. 

First  sldelobes  levels  are  shown  In  the  following  table  for  each  of  the 
first  five  PFM  waveforms,  for  two  choices  of  \|»,  namely  x/2  and  the  value  (also 
listed  In  the  table)  which  was  chosen  for  spectral  flatness.  This  table  may 
be  compared  directly  with  Table  V-1. 

TABLE  VI- I 
PFM  SIDELOBE  LEVELS 


Waveform 

SL(ii/2) 

SL(i>) 

'1' 

PFM^ 

-23.0  dB 

-18.2  dB 

5ii/8 

PFM^ 

-20.9  dB 

-17.8  dB 

19ii/32 

PFM3 

-19.2  dB 

-16.7  dB 

37.11/64 

PFM, 

4 

-18.0  dB 

-16.0  dB 

9it/16 

PFM3 

-17.3  dB 

-15.1  dB 

9x/16 

These  sldelobe  levels  are  plotted  In  Fig.  VI-15. 

The  steady  Increase  In  sldelobe  level  with  Increasing  n  represents 
another  tradeoff  In  spectral  properties,  this  time  between  near-ln  and 
asymptotic  spectral  behavior. 

The  cumulative  spectra  of  the  PFM  waveforms  show  the  same  kind  of 
variation  with  ip  as  do  those  of  the  four  waveforms  of  Section  V.'  In  Fig. 
VI-16,  the  cumulative  spectra  of  all  five  PFM  waveforms  are  shown,  all  for  the 
value  \|)-ir/2.  Figure  VI-17  Is  similar,  but  In  each  case,  the  ip-value  given  In 
Table  VI-1  Is  used,  hence  this  plot  compares  spectra  which  are  very  similar 
(and  flat)  In  the  near-ln  portion.  The  growth  of  near-ln  sldelobes  Is  clearly 
seen,  together  with  an  Increase  In  the  value  of  frequency  at  which  asymptotic 
behavior  begins  to  be  evident. 
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Fig.  VI-16.  Cumulative  spectra  of  PFM  waveforms 


It  is  interesting  to  consider  the  effect  of  letting  n  become  very  large. 

In  this  limit  the  frequency  modulation  pattern,  (^(t),  approaches  a  delta 
function  at  mid-chip,  and  the  phase  will  make  an  abrupt  jump  by  ±il).  By 
forcing  more  and  more  smoothness  on  i)i(t)  at  the  chip  boundaries,  we  have  been 
led  to  a  discontinuity  (in  the  limit)  at  t-A/2.  This  limiting  waveform  is  just 
IPSK,  with  displaced  chip  boundaries,  which  has  an  asymptotic  spectral  density 
proportional  to  f  i  The  apparent  paradox,  that  spectra  behaving  like  f 
approach  IPSK,  is  resolved  by  the  growth  of  all  sldelobes  with  increasing  n, 

and  an  increasein  the  values  of  f  where  "asymptotic  behavior"  really  begins. 

-2 

In  the  limit,  the  sldelobes  continue  high  (in  fact  going  like  f  )  for  all 
reasonable  values  of  f. 

Our  final  plotA.  in  this  section  compare  waveforms  with  equal  values  of  n. 
In. Fig.  VI-18,  cumulative  spectra  of  waveforms  A^,  and  PFM^  are  shown,  all 
for  i|)-n/2,  and  in  Fig.  VI-19,  SFSK  and  PFM^  are  compared  in  this  way,  also  for 
||)••n/2.  It  is  clear  that  the  differences  in  spectra  of  waveforms  having  the 
same  values  of  both  n  and  t{;  are  small,  and  this  is  also  true  when  these 
waveforms  are  compared  with  shared  values  of  other  than  tt/2.  These  results 
appear  to  justify  our  assumption  that  the  PFM  waveforms  provide  a 
rei>resentatlve  set  with  which  to  explore  the  spectral  properties  of  binary  FM 
waveforms  in  general. 
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VII.  WAVEFORM  COMPARISONS 


If  a  given  waveform  has  been  selected  for  a  particular  application,  the 
final  choice  which  must  be  made  concerns  chipping  rate.  It  Is  a  mistake  to 
choose  the  chipping  rate  equal  to  the  nominal  available  bandwidth,  since  this 
may  not  optimize  performance.  Against  an  optimizing  noise  jammer,  the 
performance  of  a  waveform  (coupled  with  a  matched-filter  receiver)  Is  measured 
by  its  noise-equivalent-bandwldth,  or  NEB.  If  the  peak  of  the  spectral  density 
of  a  particular  spreading  modulation  occurs  at  frequency  fj^,  then  by  the 
definition  already  given, 

G(f  )  •  NEB  -  /  G(f)  dw/2ir  . 

m  ' 


By  our  normalization,  the  right  side  Is  just  unity,  and  also 

G(f  )  -  Ag(0^)  , 
m  m 

where  6  Is  thr  corresponding  argument  of  the  spectral  peak.  Therefore  we  have 
in 

g(e„)  A  NEB  -  1  , 
or 


NEB 

^chlp 


1 

n 


For  many  of  the  waveforms  discussed  In  this  study,  6^  occurs  at  the  origin, 
and  In  any  case  this  ratio  Is  an  Important  attribute  of  the  waveform.  This 
quantity,  which  we  call  the  "chipping  factor",  cisasures  the  amount  of  NEB  (and 
hence  performance)  and  we  get  per  megahertz  of  chipping  rate.  In  this 
discussion,  system  bandwidth  Is  unrestricted. 

For  BPSK,  the  ratio  of  NEB  to  chipping  rate  Is  unity,  and  no  waveform 
discussed  In  this  study  exceeds  the  value  1.3.  In  all  the  waveforms  studied 
here,  values  greater  than  or  equal  to  unity  can  be  attained,  by  choosing  a 
value  of  i|»  which  yields  a  flat  spectral  density  In  the  malnlobe.  For  FSK,  this 
is  Illustrated  In  Fig.  VII-1,  In  which  NEB/f^j^^^p  Is  shown  as  a  function  of  \|). 

It  should  be  noted  that  beyond  a  certain  value  of  ip  (namely  1.9905  radians). 
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for  FSK 


the  spectral  peak  moves  away  from  the  origin.  For  FSK,  the  peak  value  attained 
Is  exactly  unity.  Similar  results  are  found  for  all  the  waveforms  and  Fig. 
VlI-2  shows  the  same  data  for  several,  Including  IPSK,  which  yields  the  highest 
ratios  found. 


To  choose  ®  given  waveform,  some  spectral  constraint  must  be 


Imposed,  otherwise  NEB  can  be  made  Indefinitely  large.  The  chipping  rate  Is 
simply  a  scale  factor  on  the  frequency  axis,  since 


0 


^  “  ^chlp  2ir  ’ 


and  the  spectra  depend  only  on  6.  One  way  of  recognizing  a  spectral  constraint 
Is  embodied  In  our  spreading  efficiency  parameter,  discussed  In  Section  VI.  In 
that  measure.  It  Is  simply  assumed  that  all  energy  outside  a  fixed  band  Is 
filtered  out.  More  realistic  models  should  be  compared  to  the  waveform  spectra 
In  each  Individual  case,  before  a  decision  Is  made  for  a  particular  system 
application. 

When  two  different  spreading  modulations  are  being  compared.  It  makes 
sense  to  compare  them  with  chipping  rates  chosen  to  equalize  the  NEB's,  rather 
than  with  equal  chipping  rates.  In  this  way  one  Is  comparing  the  spectra  of 
two  or  more  waveforms  having  essentially  equal  performance,  and  It  can  then  be 
judged  which  better  fits  the  spectral  constraints  of  the  problem. 

Equalizing  the  NEB  values  tends  to  equalize  the  shapes  of  the  malnlobe 
portions  of  the  spectra  of  different  waveforms  for  a  given  value  of  For 
example.  In  Fig.  VII-S,  the  spectra  of  the  first  five  PFM  waveforms  are 
compared,  all  with  i{)-Tr/2,  and  all  with  the  same  chipping  rates,  arbitrarily 
taken  to  be  100  MHz.  Spectral  density  Is  plotted,  In  dB,  relative  to  Its  value 
at  f>0.  The  actual  spectral  densities  at  the  origin  are  different,  of  course, 
since  each  waveform  has  a  different  NEB.  The  spectra  of  these  same  five 
waveforms  (still  with  \|»-Tr/2)  are  shown  again  In  Fig.  VII-4,  but  with  chipping 
rates  adjusted  to  equalize  the  NEB  values.  For  PFM]^  (MSK  In  this  case),  the 
chipping  rate  has  been  kept  equal  to  100  MHz,  and  the  others  adjusted 
accordingly,  as  llted  on  the  figure,  to  attain  the  common  NEB  value  of  61.685 
MHz.  These  waveforms  now  have  equivalent  performance  (If  they  are  not  band 


f(MHz) 


Ing  rates 
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Fig.  VII-4.  Relative  spectral  densities,  equal  NEB 


limited),  but  with  the  same  differences  already  noted  concerning  sldelobe 
levels  and  asymptotic  spectra.  However,  the  figure  shows  another  method  of 
comparison,  namely  chipping  rate  for  a  given  NEB.  The  same  kind  of  comparison 
is  made  in  Fig.  VII-5  (equal  chipping  rates)  and  Fig.  VII-6  (equal  NEB)  for 
the  five  PFM  waveforms,  this  time  with  tp-values  chosen  from  Table  VI-1,  which 
identify  "flat"  spectra  in  each  case. 

Cumulative  spectra  for  the  five  PFM  waveforms  are  given,  on  an  equal  NEB 
basis,  in  Fig.  VII-7  (il»*7r/2  in  all  cases)  and  Fig.  VII-8  (Table  VI-1 
t|;-values).  The  corresponding  chipping  rates  are  the  same  as  those  shown  on 
Figs.  VII-4  and  VII-6.  The  effect  of  equalizing  the  NEB  values  on  malnlobe 
shape  is  seen  again  in  these  figures,  which  are  analogous  to  the 
equal -chlpplng-rate  cumulative  spectra  of  Figs.  VI-16  (i|»“ir/2)  and  VI-17 
(Table  VI-1  \|)-value8).  All  of  these  cumulative  plots  show  an  inversion  in 
ordering  among  the  PFM  waveforms,  in  that  PFM^  has  the  lowest  first  sldelobe, 
and  the  slowest  asymptotic  decay,  while  higher-order  PFM  waveforms  have 
progressively  higher  sldelobes  and  faster  asymptotic  decay. 

Significant  differences  in  malnlobe  shape  appear  when  waveforms  with 
different  i|>-values  are  compared,  even  when  they  are  equalized  with  respect  to 
NEB.  In  Fig.  VII-9,  relative  spectra  are  shown  for  three  waveforms,  all 
adjusted  to  the  same  NEB.  The  waveforms  are  BPSK,  MSK  (PFM^^  at  i|;-Tr/2)  and 
"flat  FSK"  (PFMj^  at  i|;-5ti/8).  The  NEB  is  100  MHz  in  all  cases,  achieved  with 
fchip“100  MHz  for  BPSK,  and  with  the  chipping  rates  shown  in  Fig.  VII-10  for 
the  other  waveforms.  Note  that  "flat  FSK”  attains  the  same  NEB  performance  as 
unflltered  BPSK  with  essentially  the  same  chipping  rate,  but  with  superior 
spectral  shape  out  of  the  central  band  of  frequencies.  MSK,  on  the  other 
hand,  requires  a  higher  chipping  rate  to  match  the  NEB  value,  and  exhibits  a 
spectral  decay  less  rapid  than  "flat  FSK",  especially  before  the  first 
sldelobe.  The  same  comparison  is  shown  in  Fig.  VII-10  in  terms  of  cumulative 
spectra,  which  emphasizes  the  differences  in  spectral  behavior  outside  the 
malnlobe. 

This  kind  of  comparison  is  continued  in  Figs.  VII-11  through  VII-18  for 
the  other  PFM  waveforms.  In  each  case  we  compare  BPSK,  PFM^  ('|>*ir/2)  and  PFM^ 
(Table  VI-1  tj»-value8),  in  both  relative  spectra  and  cumulative  spectral  form. 
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Fig.  VII-5.  Relative  spectral  densities,  equal  chipping  rates 


Fig.  VII-6.  Relative  spectral  densities,  equal  NEB 
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VII-8.  Cuinulntlve  PFM  spectra,  equal  NEB 


Fig.  VII-9.  BPSK  and  PFM^  relative  spectra,  equal  NEB 
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Fig.  VII-11.  BPSK  and  PFM_  relative  spectra,  equal  NEB 
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Fig.  VII-16.  BPSK  and  PFM,  cumulative  spectra,  equal  NEB 
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Fig.  VII-17.  BPSK  and  PFM,.  relative  spectra,  equal  NEB 


NEB  and  BPSK  chipping  rate  la  100  MHz  In  all  these  cases,  and  the 
corresponding  PFM  chipping  rates  are  shown  on  the  figures  Illustrating 
cumulative  spectra.  The  following  table  summarizes  the  different  values  of 
NEB/f  .  .  which  govern  the  chipping  rates  used  on  these  curves. 

Cnlp 

TABLE  VII-1 

NEB/f  .  .  Values 
chip 


Waveform 

iji-value 

^’"“/^chlp 

BPSK 

- 

1 

PFM^ 

ir/2 

0.61685 

PFM^ 

it/ 2 

0.68269 

PFM^ 

ir/2 

0.72203 

PFM, 

4 

if/2 

0.74901 

PFMj 

ir/2 

0.76905 

PFMj 

Sir/ 8 

0.96353 

PFM^ 

19ir/32 

1.00847 

PFM^ 

37Tr/64 

1.02488 

PFM, 

4 

9Tr/16 

1.00502 

PFMj 

9Tr/16 

1.04070 

In  Table  VII-1  the  NEB/f(.hlp  values  for  the  PFM  waveforms  with  Table  VI-1 
Ip-values  do  not  make  a  smooth  sequence,  because  these  ip-values  were  essentially 
eyeball  choices,  and  values  of  NEB* A  near  unity  are  attained  In  all  cases. 
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VIII.  SUMMARY  AND  CONCLUSIONS 

The  motivation  for  the  study  of  spreading  modulations  was  discussed  In 
some  detail  In  Section  I.  The  chief  objective  Is  to  utilize  an  Instantaneous 
bandwidth  which  Is  much  greater  than  the  Information  rate,  and  Is  thus 
diametrically  opposed  to  the  objective  of  bandwidth  conservation  In 
conventional  data  transmission  systems.  Spreading  Is  achieved  by  subdividing 
Information  bits  Into  many  chips,  coded  with  dummy  bits  which  are  known  (In 
effect)  to  the  receiver,  but  appear  random  to  a  potential  jammer. 

Performance  measures  for  spread-spectrum  systems  are  different  from  those 
of  ordinary  data  communications  systems,  although  constant -envelope  waveforms 
with  compact  spectra  are  desirable  In  both  cases.  Sapid  spectral  decay  outside 
the  principal  band  Is  advantageous,  and.  In  the  case  of  spreading  modulations, 
spectral  flatness  within  this  band  leads  to  high  performance.  The  spread- 
spectrum  designer  wants  high  nolse-equlvalent-bandwldth  (which  determines 
performance  against  an  optimizing  jammer)  at  the  lowest  possible  chipping  rate, 
while  the  communicator,  considering  the  same  waveform  for  data  transmission, 
wants  to  maximize  bit  rate  for  a  given  allowed  bandwidth.  These  Issues  are 
discussed  In  this  report  In  terms  of  two  measures  of  performance,  the  so-called 
spreading  efficiency  and  the  chipping  factor.  In  all  cases,  performance  Is 
related  to  processing  gain,  which  depends  strongly  on  the  spectral  density  of 
the  waveform  In  question. 

A  classification  system  Is  Introduced  for  spreading  modulations,  and 
several  familiar  examples  are  Illustrated  In  their  appropriate  place.  This 
system  makes  It  easy  to  Invent  new  waveforms,  and  a  number  of  these  are  given 
In  the  report.  A  simple  method  of  computation  of  spectra  Is  developed,  based 
upon  a  representation  theorem  which  allows  all  of  the  waveforms  studied  to  be 
described  In  terms  of  sums  of  phase-modulated  pulses,  formally  very  similar  to 
BPSK.  This  permits  the  derivation  of  spectra  of  constant -phase  waveforms 
(which  Is  not  complicated)  to  be  extended  directly  to  the  large  class  of  so- 
called  binary  FM  waveforms. 

Examples  of  binary  FM  waveforms,  like  FSK  and  SFSK,  which  have  been 
discussed  In  the  literature,  are  Introduced  In  Section  V,  and  their  spectra 
Illustrated.  A  general  relationship  between  the  spectral  properties  far  from 


100 


im 

>  ’  * 


n 


K 


band  center  and  the  smoothness  of  the  frequency  modulation  pattern 
characteristic  of  the  waveform  is  derived.  This  relation  suggests  a  "graded 
sequence"  of  waveforms,  with  controlled  asymptotic  properties,  and  these  are 
Introduced  in  Section  VI.  A  "smoothness  index",  n,  characterizes  the  members 
of  this  family. 

Another  basic  parameter,  called  tj),  is  the  magnitude  of  the  total  phase 
change  per  chip  in  the  binary  FM  waveforms.  Waveforms  used  for  information 
transmission  are  usually  characterized  by  the  value  but  we  find  that 

spreading  performance  is  often  much  better  when  other  values,  intermediate 
between  7r/2  and  ir,  are  used.  Each  spreading  waveform  is  described  by  an 
Instantaneous  frequency  modulation,  and  tj;  is  basically  a  multiplicative  factor 
in  this  modulation,  hence  one  automatically  has  a  family  of  waveforms, 
parametrized  by  ip,  whenever  a  frequency  modulation  pattern  is  given.  This 
allows  generalization  of  all  the  familiar  modulations,  such  as  SFSK,  Just  as 
FSK  waveforms  are  generalizations  of  MSK.  Each  of  the  "graded  sequence" 
waveforms,  called  polynomial  FM  modulations,  is  actually  a  family,  with  ip  as 
parameter. 

It  is  shown  that  the  spectral  properties  of  spreading  modulations  of  the 
binary  FM  type  are  largely  determined  by  the  values  of  the  two  parameters,  n 
and  Ip.  The  actual  chipping  rate,  of  the  waveform  must  be  considered  a 

third  parameter,  which  acts  as  a  scale  factor  on  the  frequency  axis,  stretching 
or  compressing  the  spectrum.  When  operating  within  a  given  allotted  frequency 
band,  an  optimum  choice  of  chipping  rate  can  be  found  for  any  waveform. 

A  number  of  tradeoffs  appear  in  the  results  of  this  report.  Noise- 
equivalent  -bandwidth  versus  chipping  rate,  and  sldelobe  levels  versus 
asymptotic  spectral  behavior  are  typical  examples.  These  are  discussed  at  some 
length,  and  many  of  the  figures  have  been  designed  to  Illustrate  them. 

Hardware  Issues  are  not  discussed  here,  although  large  variations  In 
implementation  difficulty  exist  among  these  waveforms.  The  reasoning  Is  simply 
that  the  spectral  properties  of  a  waveform  are  fixed,  while  the  difficulty  of 
using  It  continually  changes  with  development  of  the  requisite  technology. 
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APPENDIX  A 


MEASURES  OF  SPREADING  PERFORMANCE 


We  consider  a  given  pulse  (or  signal  segment)  which  Is  represented  by  a 
transmitted  waveform  modulation  s(t,  b),  which  Is  T  seconds  In  duration.  This 
waveform  Is  characterized  by  the  bit-sequence  b  ■  bj  b2  ...  bj>j  (l.e.,  a 
spreading  code)  and  may,  for  example,  be  a  constant-envelope  waveform  such  as 
BPSK  or  FSK.  The  function  s(t,  b)  represents  the  complex  modulation  of  some 
carrier,  and  we  assume  It  Is  normalized  as  follows.  Independent  of  the 
bit-pattern: 

T 

/  ls(t,  b)l2  dt  -  1 

o 

The  spreading  code  changes  frequently,  and  Is  always  known  to  the 
receiver,  but  we  constrain  the  jammer  to  a  knowledge  only  of  the  spectral 
density  of  the  family  of  waveforms  (ensemble  average  over  random 
bit-sequences),  as  might  be  obtained  by  averaging  many  observed  spectra.  The 
jamming  waveform  Is  modelled  as  stationary  random  noise  with  unlimited 
spectral  agility  (l.e.,  the  jammer  Is  capable  of  generating  noise  with 
arbitrary  power  spectral  density),  but  finite  total  power. 

The  received  waveform  Is  modelled  as 

Z(t,  b)  -  A  el«  s(t,  b)  +  N(t) 
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where  N(t)  is  complex  noise  with  spectral  density  0^(1),  i.e,, 


E  N*(t)  N(t')  -  /  Gn(£)  e2»if(t’-t)  df  . 

««0 

The  total  noise  power  Is  Pj]: 

€S 

which  Includes  jamming  and  receiver  noise  (also  assumed  to  be  finite),  and  we 
define  the  normalized  noise  spectral  density 

8N(f)  -  GN(f)/PN  » 

which  obviously  integrates  to  unity* 

The  average  power  of  the  received  signal  Is 

Ps  -  a2/T 

(averaged  over  time,  not  ims),.  and  the  Input  SNR  may  be  taken  to  be 

Ps  a2 

(SNR)in  -  —  -  - 

Pn  PnT 

We  assume  that  the  receiver  Is  synchronized,  but  does  not  know  carrier 
phase.  It  Is  modelled  as  a  filter  whose  output  at  the  time  of  decision  Is 

z(b)  -  /  h*(t,  b)  Z(t,  b)  dt  . 


lOA 


The  Impulse  response  of  this  filter  depends  on  the  known  blt-pettem,  and 
decision  is  based  on  the  magnitude  of  s.  This  filter  output  contains  signal 
and  noise: 


s(b)  •  s(b)  +  n(b) 


where 

m 

s(b)  -  Aei«  /  h*(t,  b)  s(t,  b)  dt 


and 

to 

n(b)  -  /  h*(t,  b)  N(t)  dt  . 

The  system  output  SNR  is  defined  in  an  analogous  way: 


(SNR)out 


lE  8(b)|2 
E|n(b)|2 


where  the  ensemble  average  symbolized  by  E  includes  the  random  bits  as  well  as 
the  noise  N(t). 

In  terms  of  the  Fourier  transforms 


S(f.  b)  -  /  s(t,  b)  e-2irlft 

and 

00 

H(f,  b)  -  /  h(t,  b)  e-2irlft  jt  , 
•>00 

we  have 
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and 


E  8(b)  -  A  b)  S(f,  b)}  dt 


E  |n(b)|2  -  //  E{h(t,  b)  h*(t',  b)}  E{N*(t)  N(t')}  dtdt' 


Pn  /  8N(f)  E|H(f,  b)|2  df. 


Therefore, 


(SKR)out  ■  —  • 


/  E{H*(f.  b)  S(f.  b))  df|2 


/  gN(f)  E|H(£,  b)|2  df 


To  obtain  our  first  performance  parameter  we  now  assume  that  the  receiver 
uses  a  whlte-nolse  filter  matched  to  the  transmitted  waveform,  but  truncated 
outside  a  symmetrical  band  of  width  B.  This  truncation  Is  our  means  of 
recognizing  a  spectral  containment  constraint,  and  It  Is  Immaterial  whether 
the  band  limiting  Is  effected  In  the  transmitter  or  the  receiver.  We  retain 
our  definition  of  the  Input  SNR,  which  uses  the  unlimited  signal  waveform,  so 
that  the  processing  gain,  defined  as 


(SNR)  out 
(SNR)  In 


reflects  the  loss  of  signal  energy  due  to  bandllmltlng. 


The  receiver  filter  Is  thus  given  by 

S(f,  b)  ;  Ifl  <  B/2 

H(f,  b)  - 

0  ;  Ifl  >  B/2 

and  therefore 

E{H*(f,  b)  S(f,  b)}  -  E|H(f,  b)|2 

g.(f)  ;  Ifl  <  B/2 
0  ;  Ifl  >  B/2 

where 

ggCf)  =  E|S(f,  b)|2 

Is  the  signal  power  spectral  density  (ensemble  average  over  bit  sequences). 
This  density  Is  normalized  to  unity  over  all  frequencies: 

m  m 

/  88<f)df  -  E  /  |s(t,  b)|2  dt  -  1  , 

by  our  original  convention  for  the  signal  waveform. 

The  processing  gain,  using  this  (whlte-nolse)  matched  filter  against  the 
noise  spectral  density  g^^Cf),  Is  now  found  to  be 

PG  ■  T  • 

I  8N<f)88(f)  df 

B 


where  the  Integrals  are  to  be  carried  out  over  a  symmetrical  band  of 
frequencies  of  width  B. 


The  janaer  Is  now  allowed  to  optimize  his  spectral  density,  keeping  the 
total  jamming  power  fixed.  His  best  choice  will  maximize  the  denominator  In 
our  expression  for  PG,  which  occurs  when  he  concentrates  all  his  power  at  the 
frequency,  fjQ  at  which  gs(f)  achieves  the  peak  of  Its  values  within  the  band 
B.  This  optimum  choice  provides  a  lower  bound  to  the  processing  gain  obtained 
with  the  given  signal  class  and  receiver  design: 

PG  >  T 

Max  ga(f) 

B 


Now  the  nolse-equlvalent-bandwldth  of  the  band-limited  signal  Is  defined 


by  the  relation 

fga(0  df 

B 

(NEB)b  =  - 

Max  ggCf) 

B 


The  corresponding  nolse-equlvalent-bandwldth  of  the  original  waveform  would 


be 


/  8s(f)df 

m 

(NEB),  -  - 

Max  gsCf) 
f 


1 

Max  g8(f) 
f 


In  terms  of  (NEB) 3,  we  have 

PG  >  T(NEB)b  •  /  g8(f)  df. 
B 
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Because  of  our  nomallzatlon,  the  factor 
Fb  =  I  gsCf)  df 

is  Just  the  fraction  of  the  energy  of  the  original  waveform  which  lies  within 
the  band  B.  Normalizing  our  processing  gain  bound  to  the  nominal  value  TB,  we 
have 


PG  >  HqTB 

where 

(NEB)b 

ho  -  -  Fg  <  1  . 

B 

This  parameter  provides  the  desired  measure  of  waveform  performance,  given 
that  the  receiver  uses  a  white^noise  matched  filter. 

Another  parameter  is  obtained  If  we  allow  more  freedom  in  the  receiver 
design.  Specifically,  we  let 


H(f,  b)  - 


C(f)  S(f,  b)  ;  Ifl  <  B/2 

0  ;  Ifl  >  B/2  , 


where  C(f)  is  a  real  positive  function  of  frequency  to  be  chosen  presently. 
This  filter  may  be  thought  of  as  a  truncated  version  of  a  colored-noise  filter 
matched  to  S(f,  b),  where  the  noise  spectral  density  assumed  Is  proportional 
to  lC(f)]-l. 
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The  processing  gain  for  a  given  jamaer  Is  now  given  by 


I  /  C(f)  g8(f)  df  I  ^ 

PG  -  T  .  - - - 

/  C2(f)  gN(f)  g8(f)  df 
B 

Again  the  optimum  Jamming  strategy  Is  the  use  of  a  line  spectrum  at  f^,  the 
frequency  at  which  the  product  c2(f)  g8(f)  Is  a  maximum,  within  the  band  B. 
The  corresponding  lower  bound  to  the  processing  gain  Is 


PG 


>  T  .  - ^ 


/  C(f)  g8(f)  df 
B 


C2(f„)  g8(f„) 


C(f)  /  g8(f) 

/  -  —  •  /ITTO  df 

B  C(f„)  /ggCfm) 


I  2 


The  factor 

C(f)  ✓  gsCf) 

C(fm)  /g^) 


Is  an  arbitrary  real  positive  function  within  the  band,  except  that  It  Is 
bounded  from  above  by  unity.  By  choosing  this  function  to  be  Identically 
unity,  we  obtain  the  filter  for  which  the  processing  gain  lower  bound  Is 
maximized.  This  "mlmlnax”  filter  Is 
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H(f,  b)  - 


[g8Cf)]"V2  s(f.  b)  ;  |f|  <  B/2 
o  ;  |f|  >  B/2 


and  the  bound  obtained  la 


PG  >  T 


/  ^  88(f)  df 

B 


=  n  TB  . 


This  corresponding  efficiency  parameter  Is  now 

K 


h  -  - 
B 


/  /mo  df 

B 


By  the  Schwarz  Inequality,  r)  Is  never  greater  than  unity,  since 


/  /m^)  df 

B 


<  /  g8(f)df  .  /  df  <  B  , 
B  B 


In  view  of  the  normalization  of  ggCf).  It  Is  also  clear  that  n  would  achieve 
this  bound  for  a  signal  spectrum  which  Is  flat  within  B  and  zero  elsewhere, 
which  Is  of  course  Incompatible  with  the  assumption  that  8(t,  b)  vanishes 
outside  the  Interval  [0,  T] . 

It  Is  Interesting  to  note  that  the  mlnlmax  filter  has  a  flat  amplitude 
characteristic,  on  the  average,  and  hence  the  noise  power  received  Is 
Independent  of  the  noise  distribution  within  the  band  B.  The  noise  spectral 
density  can.  In  fact,  be  flat: 


gN(0 


1/B  ;  |f|  <  B/2 
o  ;  |f|  >  B/2  , 
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without  changing  the  processing  gain,  which  we  may  call  PG]^: 


PGi  =  T  /  df 

B 


If  we  knew  that  the  noise  would  be  flat,  we  would  use  the  whlte-nolse  band- 


limited  matched  filter,  say  H2: 


H2(f,  b)  H 


S(f,  b)  ;  |f|  <  B/2 
0  .  ;  |f|  >  B/2 


and  obtain  the  higher  processing  gain  PG2: 


PG2  -  TB  /  gsCf)  df 
;  B 


according  to  the  formula  derived  earlier,  PG,  cannot  exceed  PG2  by  the 
Schwarz  Inequality,  as  Just  demonstrated  above. 


However,  If  we  use  H2  and  then  let  the  jammer  optimize  against  us,  he 


will  pick 


gN(f)  -  «(f  -  fo)  , 


where  fj,  Is  at  the  peak  of  gg(f)  within  B,  and  the  processing  gain  will  be 
reduced  to  PG3: 


PG3  -  T 


Max  ggCf) 
B 


/  8s(f)  df 

B 


'  •  •  -  *  .  •  -  *,*• 


as  derived  before.  Since  the  mlnlmax/f liter  optimizes  against  an  optimizing 
Jammer,  PG3  must  be  less  than  or  equal  to  PG^,  and  this  follows  from  the 
Inequality 


/  g8(f)  df  <  Max  /  gfiCf)  .  /  /  gsCf) if  . 
BBS 

With  the  mlnlmax  filter  the  jammer  can  do  no  better  than  with  white 
noise,  but  the  receiver  does  not  use  the  whlte-nolse  matched  filter.  All 
these  distinctions  disappear  for  the  “Ideal  signal",  whose  spectral  density  Is 
flat  within  the  allotted  band  and  negligible  elsewhere.  In  this  limit, 

PGj  “  PG2  “  PG3  “  TB. 

When  we  compare  spreading  modulations,  the  chipping  rate  Is  a  parameter 
for  each  case,  hence  It  Is  natural  to  consider  efficiency  as  a  function  of 
chip  rate  for  each  waveform,  such  as  BPSK,  MSK  etc.  This  comparison  Is 
facilitated  by  Introducing  the  dimensionless  variable,  6  =  2irfA,  In  place  of 
frequency,  where  A  Is  the  chip  duration  (l.e.,  ■  fj.  Is  the  chipping  rate), 

and  the  dimensionless  parameter 


X  =  BA  -  B/fc  . 

We  define 

1 

g(0)  =  -  ggCf)  , 

A 


so  that 
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'■I  tjii  lii ; 


%  .N 


1 


—  /  g(e)  de  -  1. 

itt 


This  quantity,  g(6),  Is  the  same  as  the  normalized  spectral  density  used  In 
the  main  text.  The  matched-filter  efficiency  parameter,  no>  Is  expressed  In 
terms  of  g(6)  by  the  formula 


X  Max  g(e) 
I  e  I  <  irx 


In  the  last  step  we  have  assumed  that  gs(f)  Is  symmetrical.  For  each  type  of 
spreading  modulation,  g(6)  Is  a  unique  function  of  6,  and  the  effect  of 
variation  of  chip  rate  for  a  given  bandwidth,  B,  Is  represented  by  the 
dependence  of  tiq  on  the  dimensionless  parameter  x. 


IIA 


The  mlnlmax  efficiency  paraneter,  n»  Is  also  expressible  In  the  simple 


form 

1  • 

Figure  A. la  provides  a  comparison  of  BPSK,  MSK  and  FSK  for  the  special 
value  ip  -  1.9905t  using  the  matched  filter  measure  of  efficiency.  Figure  A. lb 
Is  a  similar  comparison  using  the  mlnlmax  filter.  Note  that  the  mlnlmax  filter 
provides  more  processing  gain  against  the  optimizing  Jammer  In  all  cases,  but 
that  the  difference  Is  slight  for  the  special  FSK.  The  difference  would 
vanish  altogether  for  an  Ideal  spectrum,  flat  within  a  band  and  zero 
elsewhere.  In  Fig.  A. 2a  FSK  waveforms  are  compared  with  ip  as  a  parameter, 
using  the  matched  filter  efficiency  criterion.  Figure  A. 2b  is  a  similar  plot 
using  the  mlnlmax  filter  efficiency.  The  optimum  efficiency  occurs  for  the 
values  Ip  >  2.2  (mlnlmax)  and  ip  -  2.0  (matched  filter),  but  the  maxima  are' 


A  /  irx  _ 

n(x)  -  —  )  /  ✓  g(e) de 

X  I  o  — 

N  A  — 


broad  In  both  cases 


Fig.  A-2b.  Minimax  filter;  FSK 
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